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PREFACE 



Throughout all history, men have been curious about the ordinary 
integers, or, more preci^ly, the natural numbers 1, 2, 3, 4, • • • . Even 
little children are fascinated by numbers, as evidenced by their games 
and their verses: '*One, two, buckle my shoe, etc!' 

The story of numl^r theory begins with the ancient Greeks, for whom 
arithmetikd was the science of numbers rather than the art of computa- 
tion* (Reckoning was called logistica and was generally deemed beneath 
the dignity of mathematicians and philosophers.) Thus Euclid, the Pytha- 
goreans, and other Greek writers were familiar with prime numbers, 
perfect numbers, amicable numt^rs, and figurate numbers. 

The classical theory dealt only with the natural numbet^. The mod- 
ern theory of numbers, however, studies the properties of the system of 
rational integers 0, ±: 1 , ±2, • • • . An integer h is said to divide an integer 
a if there exists an integer k such that a — hk. We can say that a is divis- 
ible by b, or that 6 is a factor of a, also that a is a multrple of and that k 
is the quotient of a by provided that /; 0. 

A unit is an integer that divides every integer; e.g., +1 and —LA 

Crime number, or a prime, is an integer, not a unit, that is divisible only 
y itself and the units. For example, 2, 3, 5, 19, and 37 are primes. A 
composite number, or a composite, is an integer that is not zero, not a 
unit, and not a prime; for example, 4, 21, 91, and 1 1 1 are compo" 'f!S. 
All integers take the form 2n or 2n 4^ 1. An even integer is one th.r is 
a multiple of 2; any integer that is not even is odd. 

The theory of numbers differs somewhat from other fields of mathe- 
matics in several respects. In the first place, the beginner in number 
theory needs but littfe other preparatory mathematical knowledge as a 
background — the basic principles of al^bra virtually suffice. In the 
second place, despite appearing to be relatively independent of other 
fields of mathematics and despite the innocent simplicity of the state- 
ment of many theorems, numt^r theory is noted for the difficulty of its 
problems and proofs, which require considerable mathematical insight 
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and ingenuity. To illustrate the apparent simplicity of some problems, 
consider the following theorems: 

'1) Every positive integer is a sum of four squares, and fewer than 
four squares will not suffice. 

For example: 7 = 2* + P + P+ P 

22 = 4^ + 2='+ P+ P 

(2) Every integer n can be expressed in the form n = + ~ z- 
For example 7 = 2' -|- 2' — 1' 
22 = 5'-f P -2=^ 

Such relations and properties are very simple to state and easy to under- 
stand, and are even readily illustrated by specific examples, but to give 
a general prcxjf is often exceedingly difficult. 

On the other hand, some properties are rather easy to prove. For 
example: 

(1) The product of any two consecutive integers is divisible by 2. 

(2) The sum of any integer and its square is an even number. 

(3) The product of any three consecutive integers is divisible by 3. 

Try to discover a proof for each of these by yourself! 

This collection of essays explores, in a very elementary %vay, only two 
aspects of number theory, namely the primes and the f^rfect numbers. 

—William L. Schaaf 
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FOHEwnnn 



The theory of prime nunibers can Ik? a fascinating subject. One of the 
notable acliievements of Greek niaiheniaties is Euclid's proof that the 
nuniUr of primes is infinite. His proof is simple and. in the language 
of the mathematician, "elegant!" 

Another significant cornerstone of numlier theory is the Fhtidameutal 
Theorem oj Arithmetic, which states that, disregarding the order of the 
factors, a com{M)site number can be factored int<» primes in one and only 
one way. For example: 

(1) 66- 2- II 

(2) *J6 ^ 2 - 3 

(3) 1823 = .1.7- 

It should be noted that in the field of number theory there are many 
unsolved problems as well as some "theorems " which are believed to be 
tr..c, but for wliich no pnKil has as yet been given. Thus, although it 
has lieen proved that the number of primes is infinite, the problem of 
finding the next prime after any given prime still remains unsolved. 
Nor has anyone succeeded in developing a general method for finding 
even one prime number greater than a given one. Again, it i.s suspected 
that every ixisitive even integer can be represented as the difference of 
two {Kisitive primes in infinitely many ways, but this has never been 
proved. For example: 

(1) 6 = II - .5 - l:i ~ 7 = 17 n = 19 - iS r:. etc. 

(2) 8-11 - .1 ^ n b ^ h> - 1 1 1- .11 - 2.i ^ etc. 

Perhaps one of the most celebrated unsolved problems of number 
theory is the well-known Goldbach's conjecture — every even integer 
greater than 2 can In: represented as the sum of two {Kisitive primes. 
For example: 

(1) 4^^2 + 2 

(2) 10 3 -f 7 = .5 4- 5 

(3) 80=r 11 -f 19 = 7 4- 

(4) 1(10 - 1 1 f 89 - 17 4 8.S = 41 f 59 

Here is indeed a fascinating and challenging topic! 
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by Ernst Meisst-I. who succeeded in showing that the number of primes 
below If) IS 5,761.455. The DanUh mathematician Bcrtelsen continued 
these computations and announced, in I89S. that the number of primes 
b^low 10" IS 50.847.478. This represents our most extended knowledRC 
along these unes. * 

No practicable procedure is yet known for testing large numbers for 
primality. and the effort spent on testing certain special numbers has 
been enormous. For more than 75 years the largest number actually 
verified as a prime was the S9-digit number 

2-- 1 = 170.141. 183.460.469.2S1.731.687.30S.7I5,884.I05.727. 
given by the French mathematician Anatole Lucas in 1876. In 1952 
the KDSAC machine, in Cambridge. England, established the primality 
ol the much larger (79-digit) number 

180(2'"- 1)'+ 1. 

and in the same year the SWAC digital computer, in the United States 
established the primality of the enormous numbers 2"'- 1 2*"— I and 
2' - 1. the last of which is a .S86-digit number. 

A dream of number theorists is the finding of a function f(n) which 
tvill yield prime numbers for all positive integral n. Thus f(n) ^n' - 
n + 41 yields primes for all such n < 41. but /(41) = (41)'. a composite 
number. The quadratic polynomial f(n) =n»- 79n + 1601 yields pVimes 
i»r all n < 80. Polynomial functions can be obtained which will suc- 
cessively yield as many primes as desired, but no such function can be 
lijund which will always yield primes. It was about 1640 that the great 
iii.inlK?r theorist, Pierre de Fermat. conjectured that /fn) = 2=^" + 1 is 
T'!"'' fjjr all nonnegative integral n. Forn = 0. 1 . 2. S. 4 we find f(n ) = S. 

055.H7. all prime numbers, but in I7S2 Euler proved the 
conjecture false by showing that /(5) = (641) (67004 1 7). It is now gen- 
erally le t that f(n) is composite for all other values of n, although this 
has not l)ecn established. An interesting recent result along these lines 
IS the proof, by W. H. Mills in 1947. of the existence of a real number A 
«uch that the largest integer not exceeding ^^"is a prime for every posi- 
iive integer n. Nothing was shown about the actual value, nor even the 
rough magnitude, of the real number A. 

A remarkable generalization of Euclid's theorem on the infinitude of 
the primes was established by Lejcune-Dirichlet (1805-59). who suc- 
ceeded in showing that every arithmetic sequence. 

in which a and d arc relatively prime, contains an infinitude of primes. 
I he proof of this result is far from elementary. 
Perhaps ihe most amazing result yet found concerning the distribu- 
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The Prime NinnbeiV 

Howard W. Ews 

Vri}\Hmiiim 14 of li<Kik IX of Euclid s Eiemetits is csst'iuiaUy the 
c(|uivalcnt oi the itnpomtU "lundamcntal theorem ot arithmetic;* which 
%latc% that any intffrrr greater than t run, except far the order nf the 
factors, he exfnessed as a product of pritnes in one and only one umy, 
A his theorem a$!iert!i that the prime ntimbers are building bricks from 
which all other ifiteger?» may \k made. Accordingly, the prime numbers 
have received much study, and considerable effort s have bc*en sj>ent 
trying to determine the nature of their distribution in the scMpience 
of |M>sitive integers. The chief results along this line obtained in an- 
tiquity are Euclid's pnHit of the initnittide <if primes and Erau>sihenes* 
sieve int finding all primes below a given integer n, 

Euclid's priKif. in Pro{>ositton 20 of Il<K>k IX of his Elepnents, that 
the numhfr o\ prime nnmhers is infinite, has iK'en universally regarded 
by mathematicians as a model of mathematical elegance, rhe proc^f 
employs the indirect method, or redtictio ad absurduni, and rtms essen- 
tially as follows. Siip|Kise there is only a finite numlier of prime ntunbers. 
which we shall denote by a. h,--- ,k. Set P - ri, /i, • • • k. I hen P + I is 
either prime or composite. But, since a, /^ • * • . A* are all the primes, 
/* f I . which is greater than each of a, h,'" , k. cannot In? a prime. On 
the other hand, if P + \ is c(mi|K)site. it must {>e divisible by some 
prime p, fUit /; must t>e a member of the set a. of all primes, 

which means that /; is a divisor of /f&insecpieiitly. /;canmH divide P + K 
siiue /; • L I bus our initial hyp*ilhesis that the number of primes is 
finite is untenable, and the theorem is esiabliMied. 

I1ie so-caHed sieve of flrafosfhenes is a clever device noted by the 
Cireek mathematician Eratosthenes (r b.c:.) for finding all the 
pritne numlH'rs less than a given numlHrr n. One writes down, in order 
atid starting with .S. all the odd numljcrs less than n. The composite 
numfiet^s in the secpience are then sifted out by crossing off. from .1. every 
third number, ilien from the next remaining nunilu'r. every fifth 
ninnber. then from the next remaining immlK-r. 7. every seventh nuin- 
Ikt. froTU the next remaining nmn!>er, 1 1. every eleventh number and 
si> on. In the prrnens some numlKTS will lie crossed off more than cmce. 
All the remaining nmnlHTs. along with the numl>er 2. constitute the 
list of primes less than n. 

From the sieve of Eratosthenes can lie <ibtained a cumliersinne formula 
which will determine the number of primes below u when the primes 
In-low \17.ire known. I bis formtiia w;is considerably iiuprovcd in 1870 

• AiUfm^l ffom IhiwJtd K%i% An ifitHHinrtiofi the iii%iifry of Mnthnnatus (\4'H Viitk: Riiu*- 
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tion of tfie primes is the so^^llcd prime numher theorem. Suppose we 
let .'l^ denote the number of primes Mow ri. I he prime number theorem 
then says that (A^ log^ n)/n approaches 1 as n l>ecomes larger and larger. 
In other words AJn, called the density of the primes among the first n 
integers, in approximated by 1/Iog^ n, the approximation improving as 
n increases. I'his theorem was conjectured by (tauss from an examina- 
tion of a large table of primes, and was independently proved in 1896 
by the French and Belgian mathematicians J. Madamard and C. J. de la 
Vallce Poussin. 

Exten?iive factor tables are valuable in researches on prime lunnbers. 
Such A table for all tmmbers up to 24,0(10 was published by J. H. Rahn 
in ir»59. as an appendix to a book on algebra. In 1668. John Pell, of 
Fngland. extended this table up to 100,000. As a result of appeals by 
the (ierman mathematician j. H. I^mln^rt, an extensive and ill-fated 
table was computed by a Vieimc^sc schcK)lmasCer named Fclkel. Thv 
lirsf volume of FelkeFs coinputaticms. Raving factors of numbers up to 
108.000. was published in 1776 at the expense of the Austrian imjierial 
treasury. But there were very few subscribers to the volume, and so the 
treasury recalled almost the entire edition and ccmvertcd the paper into 
cartridges to lie used in a war for killing Turks! In the nineteenth cen- 
tury, ilie ccmibined efforts of Chernac, Burckhardt, Crelle, Cilaisher. 
atul the lightning calculator Dase led to a table covering all numbers 
up to 10.000.000. and published in ten volumes. The greatest achieve- 
tneiit of this sort, however, is the table calculated by J. H Kulik (1773- 
|H6:$), at the University ol Prague. Mis as yet unpublished manuscript 
is the f csnli ol a 20-year hobby, and covers all numl>ers up u$ 1 00,000.000. 
The best available factor table is that of the American mathematician 
I), N\ I-ehmer ( 1867-19S8). It is a cleverly prejiared one-volume table 
covering numliers up to 10.000.000, 

There are many unproved conjectures regarding prime nufnbers. One 
ol these is to the effec t that there are infinitely many jiairsof txvin primes, 
or primes of the form //and p f 2, like Sand 5. 1 1 and IS, and 21>and .SI. 
Another is the conjecture made by C. Ctoldbach in 1742 in a letter to 
Fuler- (fohltiach had observed that every even iiUeger. except 2. seemed 
iepreH!Ual)le as the sum of two primes. Thus 4 2 4 2, 0 ^ .M 4 .1. 
8 5 \ .S. • • . , 16 r. 13 f .1, 18 ^- n 4 7. • • • . 48 29 I 19. • , 
100 \)1 t .1. and so forth. Progress on this problem was not made until 
W)^] wlien the Russian mathematician Schnireltnann showed that every 
{Kisiti\e integer can Ih? representi*d as the sum of not more than ,^00,000 
primes! Simu-what later the Russian mathematician \'ifiogradoff showed 
that there exists a jiositive integer iV such that any integer n > iV can 
be expressed as the sum ol at most lour primes, hut the pnn)! in no way 
jKTniits us to appraise the si/e ol N, 



FOREWORD 

A dramatic device bequeathed by the ancients is the Sieve o! £nitc»- 
thenes (c. 230 B.C.), who flourished about the time of Apollonius and 
Archimedes, two of the great^t Greek mathematicians. This is a simple 
device for testing whether or not a p(»itive integer m is a prime by 
systematically striking out all composite numbers which precede m. 
How long must this process of striking out integers be continued before 
we know that m b a prime? Eratosthenes provided the answer by the 
following theorem, which thus becomes a useful test (or a prime: 
/ -M»itive integer m is a prime if it has no pc^itive prime factor 

less tfian or equal to I, where I is the greatest integer such that P is 

less than or equal to m. 

The present article not only explains the classical Sieve of Eratosthenes 
but extends the discussion to modem developments including the so- 
called random sieve, which is then related to the "prime number the- 
orem!' This theorem, anticipated by Gauss during the first half of the 
nineteenth century, was renned and improved by 1900 in the form 

lim / ^(x) \ . 

X-* ao \ x/Iog X / 

It is of interest to note that the famous Russian mathematician Tche- 
bysheff (1821-1894) succeeded in showing that for any real number 
n > 5^ thei« is always at least one prime between n and Zn — 2. 
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Mathematical Sieves 



They sift out prime numbers and similar series of integers. 
Recent research into their properties suggests that a kind 
of uncertainty principle may exist even in pure mathematics 

David Hawkins 



It is no accident that the theories of probability and statistics are 
among the most rapidly growing branches of modern mathematics. 
Science demands them. Faced with problems too complex, or too little 
understood, to solve exactly, it falls back on laws or facts that are true 
only probably, or on the average. And from physics, considered the most 
exact of sciences, we learn that at bottom nature is inescapably uncer- 
tain and chancy. 

But if we must settle for a gambler's view of the real •> orld, can we 
not console ourselves with the thought that in the abstract realm of 
mathematics certainty is always possible? As this article will indicate, 
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Figure 1 

SIEVE OF ERATOSTHENES, a small pan of which is shown here, was devised 
more than 2,000 years ago to separate prime and composite numbers. The first 
"layer" of the sieve screens out multiples of 2 from the series of integers at the top. 
Since S passes through this layer, it screens out its own multiples in the next layer. 
Numbers at the bottom are primes which have passed thmugh all previous layers: 
they will become screening numbers in their turn. No simpler methml of deriving 
primes has yet been devised. 



the answer is by no means clear. Some provinces of mathematics arc so 
difficult that, for the present at least, they must make do with rules 
which are only probably true. Even in mathematics there may be an 
uncertainty principle not utterly unlike the uncertainty principle of 
physics. 

The text erf this sermon derives not from some new and exotic kind 
of mathematics but from arithmetic. We shall discuss the classical prob- 
lem of prime numbers. These are the positive integers — 2, 3, 5, 7, II 
and SO on — which cannot be represented by multiplying two smaller 
numbers. (Numbers which can be represented by such multiplication 
— 4, 6, 8. 9, 10. 12 and so on — are called composite numbers,) Prime 
numbers have fascinated mathematicians for centuries. It was Euclid 
who proved there is an infinite number of them. Since then many bril- 
liant minds have turned to primes and have discovered a number of 
remarkable theorems concerning them. Even more remarkable is what 
has not been discovered. For example, what is the 34th prime number? 
What is the billionth? The nth? To this day there is no general formula 
to aiiswcr these questions. The only way to find the billionth prime 
would be to write down all of the first billion and take the last. As 
another example, consider the famous twin-prime problem. Pairs of 
primes such as 1 1 and 13 or 29 and 31, which are separated by only one 
other number, are known as twins. They keep turning up in the longest 
series of primes that have yet been listed. Will they continue to recur 
indefinitely? Is their number infinite? It seems probable, but no one 
has lK?en able to prove it. 

The study of prime numbers has been quite literally as much an 
experimental as a theoretical investigation. Most of the facts that have 
l>cen proved began as conjectures, based on the insfx^ction of an actual 
series of primes. Many conjectures remain, seeming more or less prob- 
ably true. Thus an indispensable tool of the number theorist is a long 
list of primes. 

One of the best known, now found in every well-equipped mathe- 
matics library, was compiled by D. N. I^hmer of the University of 
California in 1914. The volume contains a table of the 664,580 prime 
numbers smaller than 10,000,000, plus a few more to fill the last column, 
ending with the prime 10,006,721. Lehmer's work was completed before 
the age of automatic computation; today there are even longer lists, the 
longest h^ing "published" only on magnetic tape. 

Modern tables of primes are prepared by a njethod, essentially un- 
altered for 2000 years, which is called the sieve of Eratosthenes. Its 
inventor was one of those great figures of the Hellenistic Age who seem 
today, across the intervening centuries, so clairvoyant of the spirit of 
modern science. Eratosthenes of Alexandria is best known for his feat 
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In the randmn lieve. however, the specific numbers to be cUminatcd are chosen by a random Pj««» »y^*«'Jif?.*2liSlJ 
colored wheeU, Thus the random sieve produces a difercnt set of numbers each time it is used, while the set of prime 
numbers is invariant. 

Figure 8b 



of measuring the size of the earth. But he was a man of universal learn- 
ing who wrote also on geometry, the measurement of time, and the 
drama. In his own day he was nicknamed "Beta" because, it was said, 
he stood at least second in every field. Modem electronic computers can 
make far longer lists of primes than Eratosthenes could have, but his 
principle of computation has not been much improved. 

The method is almost obvious {see Figure 1 ). Simply write down a 
series of positive integers and proceed systematically to eliminate all the 
composite numbens. The numbers that remain -that fall through the 
•sieve' - are primes. We begin by knocking out the even numbers 
which are multiples of the first prime number: 2. (One is not usually 
called a prime.) When we have done this, the smallest of the remaining 
numbers is the second prime: 3. Now we eliminate the multiples of 3 
irom the numbers which survived the first sieving operation. Five is the 
next number remaining, so its multiples drop out next; then the multi- 
ples of 7, and so on. 

The reader may wish to try a somewhat longer version of the sieve 
than the one shown in the illustration, where 7 is the largest sieving 
number. In number theory the distance from the obvious to the pr<v 
found is sometimes very short, and any amateur willing to play the game 
is on the verge of some first-class mysteries. At any rate, a little manip- 
ulation of the sieve will make clear some of its properties. Every sieviriff 
number is a prime. The first number sieved out by each one is its own 
square: the first number eliminated by 2 is 4; by 3. 9; by 5. 25 and so on 
In addition, the fraction of the remaining integers eliminated by each 
sieving number is its own reciprocal: 2 sieves out half of the remaining 
numbers, 3 sieves out a third, 5 sieves out a fifth. 

By carrying out the sieving operation through the prime number 31 
we can obtain all the primes in the first 1.368 integer. (The first number 
sieved out by 37. the next prime, is 37% or 1,369.) For purposes of illus- 
tration we have arranged the first 1,024 of the integer in a 32 X 32 
array, with the prime numbers shown in italics (see Figure 4). The list 
IS short, but it does demonstrate that the frequency of primes slowly 
decreases in a rather irregular way. From considerably longer tables 
Adrien Mane Legendre. and later Karl Friedrich Gauss, were able to 
guess one of the most important facts about primes -the celebrated 
Prime Number Theorem. This tells how many primes we may expect 
to find by carrying the list out to any given number. It states that if 
the number is n, then there are about n divided by the logarithm of n 
(n/log n) primes before it. As n grows larger, the error in the formula 
becomes a smaller and smaller proportion of the exact numl^er of primes 
(iauss. whose skill in computing belied the myth that mathematicians 
cannot add and subtract, arrived at the theorem by a combination of 
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DISTRIBUTION of random "primes" between I and 1,024 (Figure random sieve might yiehl an even more similar distribution. I he 

5; resembles that of true prinies in the same number sequemt? resemblance of the two series lends to intensify a^ they are imreased 

(Figured). Both sets of numbers (in bold face) thin out irregularly as in length, 
the sequence progresses (see totab at right), Another *'run" of the 
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arithtncticat insight and purely empirical siudy. It wan nm pmved for 
almost another century. In the 1890s the Ikrlgian mathematician Charles 
de fa Vallee pQU!isifi and the French mathemattcian jactjues tiadamard 
tnde(x.*ndently foinid a pnK>f, but it made use of concepts outside simple 
whole numbers. It was not until 1950 that the Norwegian mathematician 
Allc *Sclberg discovered a purely arithmetical proof. In the quaint vocab- 
ulary of number theory his proof is called elementary, but it is nor easy, 

The difficulties of the Prime Numlnrr Theorem are connected with 
the pu//lingly irregular way in which the primes aa* distributed. Indeed, 
the theorem itself dovs no more than state a statistical average. Out- 
rageous as it may seem, the sequence of primes is just as "random" as 
many of the natural phencmiena on which we make bets. Sometimes we 
think that if we knew enough about the individual events of which such 
phenomena are compf>sed. we could predict their cmicome with cer- 
tainty. This is surely true of the primes. The sieve will e\entually tell 
us about the primality of any given numlier. Hut ir cannot tell us aboiU 
all numbers, iiecause the sequence is itself an infinite, unending process. 

From the time of (iauss mathematicians have talked. |>erha{^ rather 
shameiacedly. about the "probable" fiehavior of primes, and this kind 
of reasoning has been very helpful. No mathematician, however, seems 
to have gone itie whole way and made a purely statistical model of the 
prime numiier distribution. Receiuly I was to try it. and I found 
that the model helps clarify the Prime Numl>er Theorem. Furthermore, 
it places the whole subject in a new perspective. In particular, the the- 
orem ni* longer ap|>ears as a s|K-cial fact al>out the secjiience of nmnlKTS 
which cannot be produced by multiplying two smaller numl>ers. but 
rather as a conunon feature of all sequences of numbers generated by 
sieves of a certain ty}>e. 

The fuodel is culled the random sieve, and ii works like this (see 
Figure ^ j. .Start with 2 as the first sieving numl>er. just as in the method 
of Fratostlienes. Now make a kind of roulette wheel that is divided into 
two e(jual parts, black and white, (io down the list of integers following 
2. and for each one spin the wheel. If the black part of the wheel stops 
at the iMiinter. strike the integer out; if white stops at the |K)inter. leave 
the intej;er in. Note what you have accomplished. In the lony; rtm you 
will have sieved out half of the integers, just as the first step in the 
prime ninnber sieve does, liut jtist which ones go out is a matter of 
chance, and the list will be different each time you try it. 

Next take the first number that was not removed. Supjiose it %vas 4. 
Make a new wheel of which a fourth is black and three-fourths is white. 
Spin the wheel ffir each succeeding numlnrr left after the first sieving. 
VVhen blafk comes up. strike the mnni>er out; when white comes up, 
leave the number in. t his time you have removed a fcmrth of the re- 
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mainin^ numbers. Proceed again to the first number iiot removed — 
say 5. Repeat the procedure using a sieving probability ot i/^, and so on. 
After any number of steps you will be left with a series of integers %vhich 
might be called "random primes!' 

If you want to try the sieve yourself, you need not actually make 
roulette wheels. A ubie of random numbers or, failing that, a telephone 
btxjk. will do. Express each sieving probability as a four-digit decimal 
(r,g„ 14 = .25f)0). For each "spin of the wheel" read successive telephone 
numbers. If your probability is »4, then any number whose last four 
digits are 2499 or less tells you to eliminate the integer in question: 
2500 or more means to leave it in. 

One run of the random sieve for the first 1,024 integers is summarized 
in the table on the preceding two pages. Comparing the distribution of 
these random primes with the actual ones, we can see that our sieve 
acts something like the sieve «f Eratosthenes. This is partly in spite of 
the random element, but partly because of it. For a much longer series 
the general statistical similarity would be even closer. 

It may seem paradoxical that we can take a statistical model, involving 
an infinity of random choices, as ersatz for the straightforward and per- 
fectly defined sieve of Eratosthenes. The paradox is the same as the one 
wliith underlies statistical mechanics: the average behavior of an assem- 
bly of molecules is easier to describe than the actual behavior of any one 
i>f them. Of course the random sieve preserves only the general features 
«>f the prime fHimber sieve. The eccentricities of the latter are averaged 
out by randomizing them. In either case any number not sieved out 
iKTConies in turn a sieving number. It starts a process by which a propor- 
tion of later numbers is removed, ecjual to the reciprocal of that sieving 
number. Every wave of sieving in the prime-number sieve, except the 
first, is determined strictly by the result of previous waves. At every 
correspmding point the random sieve makes probability choices, partly 
determined by its own earlier statistical liehavior. 

How closely the random sieve actually approximates the sieve of 
Eratosthenes is demonstrated by the fact that the Prime Number The- 
orem holds for random primes. This can l)e proved by some elementary 
mathematics, which in this case is also fairly easy (srf Figure 6). 

Perhaps the parallel between the two sieves is not so surprising. We 
might say. indeed, that ihc prime-number sieve would have to be re- 
markably abnormal in its detailed behavior not to lead to the same 
general result as the random sieve. This statement implies that the ran- 
dom sieve can be taken as a criterion of normality. 

If so, there must be other sieves — in fact, an infinite number of other 
sieves — that have the same general chai^cterbtics as those of the sieve 
of Eratosthenes, but which differ somewhat in the details of their defini- 

18 



ERIC 



tion. They will not yield the prime numbers in general, but numbers 
having snnie other fifiecifil projierty. In 1956, as it hap{>ens. Siani<»taH M. 
Ulatn and his associates at the Los A!amf)s Scientific Laboratory pub 
lished some results of a new type of sieve which yields what they called 
"lucky" nunilnrrs. Their sieve begins by removing the muhiples of 2. 
leaving 3 as the first number not sieved out. Instead of removing next 
the multiples of 3, the Ulam sieve removes every third remaining muii- 
ber. Since 5 is the thitti number in the list of remaining numbers, it 
drops out, but 7 remains. Hence in the next wave every seventh number 
of those still remaining is eliminated, and so on (,see Figure S). The num- 
Inrrs that escape are "lucky;* It has been proved that the analogue of the 
Prime Numlx-r Theorem holds for lucky numbers. Thus the random 
sieve is a model for the lucky numbers as well as for the primes. 

So far the random sieve has only duplicated results that can be ob- 
tained independently and rigorously for the sieves of Eratosthenes and 
ot IJIam. The mathematics of it, however, is mostly easier. Therefore 
many additicmal theorems can lie obtained from the random sieve and 
conjectured to l>e true of the other two. Such conjectures are not proofs, 
but we can say that unless the prime number and lucky sieves are vastly 
abnormal, the results must hf)ld for them. 

Ixt us lo€>k at a couple of examples. As we go to larger and larger 
numlKTs in the table of integers, the spacing iK'twcen successive primes 
(or liickies or random primes) grows greater in an irregular way. In the 
neighlHirhood of any number. «, the average interval is about the loga- 
rithm of n. What is the greatest interval? We do not know the answer 
for primes or luckies. But for the randotn sieve we can prove that, with 
f)nly a finite number of exceptions, the interval is never greater than 
the scpiare of the logarithm of n, that is, (log n)-. Vhc chance that there 
win Ik? any further exceptions can In? made* as small as we please by 
taking a sufficiently large n. No upper boundary to the interval t^tween 
successive primes or successive luckies has been found which is anyw*here 
nearly as sfiiall. although from the existing tables it hx)ks as though the 
formula should hold for them uxy. 

Another example is the twin-prime prf)bleni mentioned earlier. In 
the randcmi sieve there is almost certainly an infinite number of twi^is. 
Indeed the average interval between twins ought to be alniul (log u)'\ 
and the fnaximmn interval iKiween them, with only a finite number 
of exceptions, ought to \w (h)g n)\ Again the tables suggest that these 
results are also true for primes and luckies, but no one has any idea how 
to prove such results. 

Although the randcmi sieve does not solve any classical probk.ns 
concerning primes, it does enable us to reformulate sue h problems. We 
may ask: "Are the prime numbers normal in "uch and such a respect?" 
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The random sieve, or certain modifications of it, defines what we mean 
by normality. If the properties we are talking about depend on the 
exact fine structure of the sequence of primes, the answer will obviously 
be no. Thus all primes except the number 2 are odd, while . this. is in- 
finitely improbable in the sequences of random primes. But average 
properties such as those we have discussed do not seem to depend on 
the fine structure, and those may be presumed to be normal for primes 
or luckies. Can anyone find a major abnormal property, in this sense, 
of the sequence of primes? Or the sequence of luckies? 

In the opinion of the author the concept of normality raises some 
very deep questions about numbers and the theory of numbers. Sieves 
as a class are a type of feedback mechanism: the output of one sUge of 
the process determines the input of the next stage. Now in any such 
mechanism the nature of the coupling between output and input is 
crucial; the result may be stable and predictable for one type of coupling 
and unstable for another. So far as the outcome of the random sieve is 
concerned, it is in one respect extremely stable. If by chance there are 
relatively few sieving numbers in the early stages, they will remove 
relatively few later on, and so there will be an increase in the later stages 
to compensate for the initial deficit. The sieves of primes and luckies 
share this characteristic. But this is a statistical stability. 

When we look at other aspects of the prime or the lucky sieve, how- 
ever, we find elements of instability. The detailed ordering of primes or 
luckies depends upon the individual sieving numbers that precede 
them, and this involves a growth of complexity without apparent limit. 
Some easily defined properties of normal sequences, for example the 
two described, may depend strongly enough on this complexity to make 
it impossible, in a finite number of steps, to prove that they hold. Here 
is the analogy, if it be one, with the uncettainty principle of physics: An 
infinite complexity requires infinite time to resolve it. If our suggestions 
have substance, we will have examples of mathematical statements which 
are almtist certain, but which cannot, in principle, be proved. Examples 
ot undecidable propositions are known in modern arithmetic [see 
"Godel's Proof;' by Ernest Nagel and James R. Newman; Scientific 
American. June. 1956], but so far none of the unproved conjectures 
about prime numbers has l^cn shown to be undecidable. Perhaps none 
of them is. If any are. however, the random sieve will be a model for the 
primes in a deeper sense than any we have exploited in this article. We 
cannot distinguish an infinitely complex order from a random one, and 
so we might be forced to admit that there is a certain background of 
noise even among the eternal verities. 
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"LUCKY*' NUMBER SIEVE resembles the prime and random sieves already described. Here, also, numbers which are 
not eliminated become sieving numbers and remove a proportion of the remaining numbers equal to their reciprocals. 
Elimination is by counting: thus S removes every third remaining numben 7 every seventh. Like primes, the lucky" 
numbers form an invariant series. 

Figure 8 



2 3 7 9 13 15 2X 25 31 ★ 33 37 43 49 51 63 # 67 ®? 73 75 79 87 «3 ♦ 99 

im III 115 . 127 it 129 133 135 141 151 159 ^ 163 169 171 159 # 193 19S 201 2C» 211 219 223 # 231 235 

237 241 #259 261 267 273 2«3 ;»5 # 259 297 303 307 319 # 321 327 331 339 349 * 357 361 367 * %5 391 

393 399 409 415 ★ 421 427 429 433 # 451 463 475 477 ^ I 467 469 496 511*517 519 529 539 537 541# 

553 959 # 977 579 563 591 601 # 613 615 619 BZl 631 639 # 643 645 691 695 # 673 679 ^3 699 * 717 

723 727 729 735 * 739 741 749 * 769 777 791 787 # 801 iOB 819 S23 631 * 841 855 # 667 873 883 889 899* 

897 903 925 n7 it m 933 937 S^7 * 961 979 979 981 991 * 993 997 10^ lOtl 1021 l<m * 




DISTRIBUTION OF "LUCKY" numbers between 1 and 1,024 resembli» that of primes and landom prim<s, thinning out 
gradually but irregularly as the list increasaL This table shows only the "luckies ; the intervening numbers are omitted. 
Stani set oflF luckier within succmivc series of $2 integers: eadi of these grou|^ corresponds to a sii^e line of the tables in 
Figures 4 and 5. 
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DERIVATION OF THE PRIME NUMBER THEOREM 
FOR THE RANDOM SIEVE 

Let us consider the fate of any two consecutive numbers, say 127 and 
128, on a run through the. randoni'Sieving operation. We shall com{^re 
their probabilities of getting through the sieve; i.e., of becoming sieving 
numbers or "random primes" themselves. 

Call these probabilities P,„ and P,^,. Now it is obvious that 128 runs 
the same risk of being eliminated by previous sieving numbers as does 
127, except for one possibility. If 127 becomes a sieving number, it can 
eliminate 128, but not vice versa. The probability that 127 is a sieving 
number is P,„. If it is a sieving number, the probability that it will elimi- 
nate 128 (or any other following number) is 1/127. The chance that the 
two events will occur and that 127 will eliminate 128 is the product of 
their probabilities: P,„ X 1/127. The probability that this will not hap- 
pen is 1 - P,„/127. Except for this factor the chance of survival for 




AVERAGE NU\fBER OF RANDOM "PRIMES" in the first n integers is shown 
bv the area under the solid curve, roughly approximated by the hatched recungle 
(drawn here for n == ISO), The area under the broken curve gives the approximate 
number of true prim^. Since the two curves approach each other as n increases, 
the two sets of primes are very like. 

Figure 7 

128 is the same as that for 127. Its net probability is therefore the product 
of the two: P,,, = P,„(I ~ P,,,/127). 

At this point it will more convenient to shift from t!ie probabilities 
to their reciprocals. The reciprocal of a probability has itself a clear sta- 
tistical meaning: it gives the average interval, or range, between two 
events. (Instead of saying that the probability of double six in dice is 



1/36, we can as well that the average interval between throws of 
double six is 36,) Denote the reciprocal of P,,, by X,^^, and of P^^, by X^^.- 
X,j,^ measures the average interval between sieving numbers in the neigh- 
borhood of 127 and X^^^ measures the same interval in the slightly shifted 
neighborhood of 128. 

By a little algebra we can show that if P,„ =: P,„(l ~ ?^^J\27), then 
= + 1/127 + r, where r is a negligibly small remainder. For 
practical purpc^es we am wiy that X,5^ = X,^^ + 1/127, Now a similar 
argument would show that X,,, = X^^. + 1/126, and so on. Eventually 
we arrive at the result that X,^ = 1 + 1/2 + 1/S + 1/4 . . . + 1/127. or. 
in general, X^ = 1 + 1/2 + 1/3 + 1/4 . . , + 1/n, with a remainder that 
is still negligibly small In calculus books we discover that the ^ries 1 + 
1/2 + 1/3 + 1 /4 . . . + 1 /n is nearly equal to log n for fairly long series. 
The difference can be made as small as we like by making n large enough. 
Therefore we can say that, in the long run, X^ = log n, or P. = 1 /log n. 

The graph on the preceding page shows the values of 1 /log n (and, for 
comparison, the reciprocal of the actual values of the series 1 + 1/2 + 
1/3 1/4 ... + 1/n). Thus the curve is also a graph of P^* Suppose we 
now want to know how many random primes, on the average, there 
should be before any number n. We simply add the probabilities that 
each smaller numtn^r Incomes a sieving number. Graphically this is the 
same as taking the area under the curve. But if n is very large, then the 
difference between the area under the curve and the area of the shaded 
rectangle, which is n X P, is negligible. Hence we can say that the aver- 
age number of random primes out to n is n X P^. But P^ = 1/log w, so 
the number becomes n/log n. And this is the Prime Number Theorem! 

Having completed the proof, we may reexamine our reasoning to see 
why the result is plausible. The essential step was to find that X^ + I = 
X^ + I In. This equation says that on the average, over many repetitions 
of the sieve, any number n removes enough of the numl^rs following to 
lengthen the interval between them by l/«. Take a specific example. Sup- 
pme that is I /5 and X^^^, is 5. Then 127 will be a sieving number 1/5 
of the time. When it is, it will eliminate about 1/127 of the remaining 
numbers, lengthening the average interval Ixrtween them from 5 to 5 + 
5/127. Since it only does this about one time out of every five trials of 
the sieve, its average effect will Ik? to letigthen the interval from 5 to 
5 + 1/127. 

The same chain of reasoning is plausible for the prime-number sieve. 
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You have already seen how the Sieve of Eratmthenes can he used to 
determine all the primes up to any desired number. In the present essay 
wc find an ingenious modification of the sieve in the form of a mechan- 
ical chart, which reveab additional properties of the primes. 

The property of greatest interest, perhaps, is the fact that any prime 
greater than 3 is equal cither to one more or one less than a multiple 
of 6. Although the author proves this property, you might like to refer 
to a table of primes and verify this property for a few cases. 
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The Factorgrain 

Kennetli P Swallow 

Thi^ problem of Kiiding all the prime numbers has intrigued mathe- 
maticians through^ the ages. The nn^y attempts to solve this problem 
have yielded only methods that will produce a finite numlx^r of primes, 
the most noted of these being the Sieve of Eratosthenes. 1 he Factorgram 
is an adaptation of this systematic mechanical method. 

In the Eratosthenes Sieve, to find all the primes less than a selected 
number, AT, all the integers from 2 to N are written in order. I'he num- 
ber 2. which is known to be a prime, is encircled and every second 
number from 2 is crossed out. llicse are the nmltiplcs of 2 and hence 
cannot lie primes. 

(D s \ 5 X 7 X 9 >i( n 

X IS H 15 >^ 17 >^ 19 ^ etc. 

I he numlKT 3, which is prime l^cause it is the only remaining number 
less than 2\ is encircled and ever)' third number from 3 is crossed out. 

(D®\5X7>&X>^11 
X H >^ H '7 X 19 ^ etc. 

Now, 5 and 7 are the only remaining numbers less than $\ therefore 
they must l>e prime numl)crs. This process is continued until every^ 
multiple of every prime number up to VA^ is crossed out. The remain- 
ing numlxTs are the prime numlnrrs less than 

HOW TO %fAKF A FAC:TC)RGRAM 

lb lind all the prime numbers less than a selected number. A', by the 
Factorgram, place all the numbers from H to iV in rows of six numt>ers 
as follows: 
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Now the nuiltiples of 2 can l>e crossed out by drawing lines through the 
entire first, third and fifth rolunms, with the exception of the numl>er 2 
itself. Similarly, the nuiltiples of 3 can \w crossed out by drawing a line 

* Krtmrth H SwalloH. "Flmtriii4i> NtimJ>ct Theory in High Sthw>l M«sitirm«iluiii' p{>. Hi-9%. I n 
piihitfthni M.iitrr'i 1 hni«, Ohio Stdtv rntvnuitVi IQfiS. 
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through thf cn/irc fourth.column. with the exception of the number 3 
itself. The first column contains muhtples of S but it is already crossed 
out. Next, the niuhiples of 5 arc to be crossed out. The first six of these, 
5. 10. 15, 20, 2.5, and 30. lie in a straight line running diagonally down- 
. \^rd from right to kft. Tht: next six multiples of five (35 to 60) lie in 
another straight line, which is parallel to the first line. All the multiples 
of 5 can be crtissed out by a set of such parallel lines. Next, the multiples 
of 7 can be crossed out by a set of such parallel lines running downward 
from left to right. The multiples of all prime numbers can be crossed 
out by similar sets of parallel lines. In the Factorgram as in the Eratos- 
thenes Sieve, when all the multiples of all the prime numbers less than 
vATare crossed out the remaining numlicrs less than N are all primes. 

The Factorgram can be made c)n a piece of paper and then rolled into 
a cylinder so that the numbers form a helical spiral. (In Figure 1, roll so 
that the two zeros coincide.) In this form, each of the sets of parallel 
lines which cross out the nnjitiples of the prime numbers will also form 
a helical spiral. 

FKATURKS OF THE FACVfORGRAM 

The main purjKise of the Factorgram, as ol' the Eratosthenes Sieve, 
is to lind all the primes up to any selected numl>cr. However, the Factor- 
gram has many features not Itnnid in the usual Sieve, 

1. File mechanical process is very easy. I he cohunns of numbers can 
be made (juickly with a ty|x'writer. If a long Factorgram is to be made, 
jieriods should Ik- placed after the numbers as was done in Figure I. 
I he i>eriod. rather than the figure, is used to represent the exact location 
of each number. (In Figure I the distniue from the zero line to each 
tiuml»er is profwrtional to the magnitude of the tiunifjer. This improves 
the Factorgram in its cylindrical form but is not really necessary for 
prcjpcr operation.) A pair of draftsman's triangles can be used to draw 
the parallel lines needed to cross out the multiples of each prime num- 
ber. The first line of each set of parallel lines is determined by zero and 
the jirime number. All such lines pass thrrnigh zero, since zero is a multi- 
ple of every number. 

2. The prime numbers, which seem to be so haphazardly scattered 
through the numlKT system have, with the exception of 2 and 3. settled 
driwn to cxc upy |>osiiions in only two ol the Factorgram s six cohunns. 

.1. The presence of prime pairs of the form p and p -f 2. such as 5 
and 7. 1 1 and 1 3, etc.. and of the fonn /; and p f 4. such as 7 rmd 11, 1.3 
and 17. etc. l>ecome more obvious. Also, the relationships of prime 
nunil>ers to the numljcr 6 are emphasized. 

\. just as a prime nmnlHT can Ix* identified by the lack of lines pass- 
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ing through ii, a composite nuiiibcr can be identified by the one or more 
lineal passing through it. These lines provide a means, tree of all trial- 
and error methods, tor Hnding all the prime factors of a composite 
number. It is this property of the Factorgram which gives it its name. 

J he method of factoring a composite number by the Factorgram is 
as follows: First locate the number on the Factorgram and trace any one 
of the lines passing through it back to its prime number origin (the 
last number on the line before reaching zero). 1 his number is one of 
the prime factors of the original comp)site number. Divide the original 
tunnl>er by this prime factor to obtain a second factor. Find this new 
factor on the Fac torgram to see whether it is prime or com{x)site. If it is 
prime, the problem is completed; if it is composite, continue the process 
until the factors of this factor are prime. 

For a numerical example, consider the factoring of 117. There are 
two lines passing through it on the Factorgram. One of these goes back 
to ^. Dividing 117 by ,S we have V.). On the Factorgram, 39 has two lines 
passing tlirough it also, one going to ^ and the other to 13. Therefore, 
the factors of 117 are 3 • 3 • 13. 

Frequently the divisions will be inmecessary because there may be 
as many lines passing through the given number as there are prime 
factors of the numiier. In this case, each line will give one of the prime 
factors. In fact, the necessity for division occurs only in two cases, 
(a) if the s;ime nuhilnrr occurs as a factor two or more times, and (b) if 
one of the factors is greater than the square root of the largest number 
cm the Factorgram. V^ariations of the Factorgram that will eliminate the 
division in both of these cases can be made but these variaticms become 
overly complicated with uyty many lines, 

fi. II two innnlKTs have a common factor, they will be ccmnectcd on 
file Factorgram by the line representing that factor. For example, 26 
and (ifi are ccmnectcd l)y the line that passes through 13; 88 and 121 
are connected by the line that passes through 11; and 70 and 105 arc 
connected by two lines, one passitig through 5 and the other passing 
through 7. This proj>erty of the Factorgram cran lyc useful in reducing 
fnic tions. By locating the numerator and the denominator of a fraction, 
one can tell whether the fraction can be reduced and if so, by what 
nuuilH-r the innnerator and dencmiinator should be divided. 

WHV IMF FAC rORGRAM WORKS 

The umisual projHTties of the Factorgram are based entirely ujKm 
the following pro|K)sition: 

llienrem: All prime mnnl)ers greater than 3 are either one more or 
one less than a tnultiple rf 6. ' 

If! fither words, all primes greater than 3 are given l>y cme of the two 
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expressions, 6n + 1 or 6n — 1. For example, 6(1) -f 1 = 7, 6(2) +1 — 13 
and 6(3) + 1-19, while 6(1) - 1 ^ 5, 6(2) -1 = 11, and 6(3) -1 = 17. 
Of course, the convene of this statement is not necessarily true. For 
many values of w, the expressions 6n + 1 and 6n — I will not yield primes; 
for example, n — 4 in 6n -f- 1 gives 25, n = 6 in 6n — 1 gives 35, and 
rt 8 in 6n -f 1 gives 49. . . - 

The proof of this theorem is c|uite simple. Every number can be ex- 
pressed by one of the following six forms, four of which are always 
factorable, if n is greater than zero. 



6n =6(n) 

6n + 1 not factorable 

6n + 2 =2(3n + l) 

Cm + 3 = 3(2n + 1) 

6n + 4 ^ 2(3n + 2) 



fm f 5 (or 6n — 1) not factorable 

I lie expressions fm + 5 and 6 (;i +1) — I are equivalent since 5 more 
than a multiple of 6 is also I less than the next multiple of 6. Obviously, 
if four of these six expressions are always factorable, the primes must 
l>e expressed by the other two expressions, and hence, the theorem is true. 

N«n\\ in the Factorgram, n is the munlK-T of each row (if the first 
row is (»), and the six columns are, from left to right, 6n, 6;i + 1, 6n 4- 2, 
6;i I 3, 6n + 4, and Gn + 5. 
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riierefore. from the above theorem, all the prime numbers above 3 
must lie in the second and fifth columns. 

The multiples of a prime of the form 6n + 1 must be 2(6n + 1)* 
3((m f 1), 4(6;i f 1), etc. These, when simplified are 12n +2, 18n + 3, 
24;i f 4, etc.. or 6;i' + 2. 6n'' + 3. 6n'" + 4, etc. Hence, the multiples 
of 6;i 4 1 must progress in regular fashion from the siTcmd column to 
the third column, the third to the fourth, the fourth to the fifth, and so 
on. In a similar manner, the multiple's of a prime of the form 6n - 1 

' In coitgt tti*sKv nofaf i«>n, for p > 5, for ff^± \ (ifuxj fi). 
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(the same as 6m -f 5) progress from one toluimi to the next, but in this 
case they progress from right to left. 

The number of rows tfie line passing through the nuihiples <if a prime 
goes dtiwn as it progresses fonvard or backward from one cohnnii to the 
next is given by the vahic of n for that prime. For example, the line 
through the multiples of 7, for which n = I in 6« + 1. goes down one 
row as it progresses fonvard one columi». while the line passing through 
the multiples of 19, for which « — 3 in 6n + 1 . goes down three rows as 
it progresses forward one colunni. If the prime is of the fonn 6« — I 
the line passing through its multiples will progress backward (right to 
left) instead of forward. 

This sort of slope is helpful both in setting up the parallel lines in 
making the Factorgram and in using the Factorgram in factoring. 
In a long Factorgram it is not necessary to trace the parallel lines or 
spirals back to the prime that produced them. One merely needs to 
locate the numl>er to lur factored and note how many rows down the 
line (or lines) gcKs as it progresses forward or backward one column. 
This number is the value of n which is to be substituted in 6m -f- 1 if it 
progresses forward or 6m — 1 if it progre.sses backward. The value of the 
resulting expression is the same prime numl>er that would be obtained 
if the line were traced back to its origin. 

While the Factorgram is neither particularly profound nor useful, it 
is simple eiiough for high school students to understand and offers many 
f>p{)ort unities for interesting classrcMMU or mathematics club discussion as 
do the Kratosthenes' Sieve. Pascal's Triangle, and Magic .Scpiares. 
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Consider" the "factors of the integer 6. namely. 1, 2, S, and 6; their sum 
is 12. or twice 6. Again, the factors of 28 are 1. 2. 4, 7, 14, and 28; their 
sum is 56. or twice 28. Such numl^rs are called perfect numbers. An 
integer is said to be a perfect number if the sum of its factors is twi( e 
the given integer. A perfect number is also defined as any integer which 
equals the sum of its proper factors, where a "proper" factor of a number 
means any of its divi^rs except the number itself. 

The ancient (irecks were familiar with perfect numbers. In fact. 
Euclid proved that if an even integer is of the form 

2p « . (2'' - 1), where 2" — 1 is a prime. 

then that integer is a jierfect number. 

The converge theorem was proved by Euler some two thousand years 
later. If an integer is an even perfect number, it has the form 
2M t . ^9P _ where both p and 2*" — 1 are primes. 

It is interesting to note that all known perfect numl^rs are even. 
Although no ixid perfect number has ever been found, mathematicians 
have not yet succeeded in proving that none exists. 

Integers of the form 2" - I. where // is a prime, are called Mersenne 
iiuinluTs. after the French niutheniatician Marin Mersenne (c. 1620). 
rhey play an important role in the study of perfect numbers. Me?senne 
nuinticrs arc designated as Mp — 2" — 1 . where is a prime. Thus, for 
p — 5. M, = 2'' ~ I = 31. If Mp is a prime numt^r. it is called a Mer- 
senne prime. Until recently, only 20 Mersenne primes were known. In 
HI6.1 the three largest known Mersenne primes, M^„,«. M,^,,. and M,,,,,,. 
were discovered by the electronic computer Illiac III at the Digital 
Computer Laboratory of the University of Illinois. This bring to 23 
the numl>cr of known Mersenne primes which are 

2. ^, .5. 7. 13. 17. 19. 31. 61. 89. H>7, 127. 521, 607. 12V9, 2203. 2281, 
3217. 42.53. 4423. 9689. 9941, 1 1.21.i. 

A multiply perfect number is an integer n the sum of whose factors 
is a multiple of n. For example, the sum of the "proper " factors of 120 
is twice 120. or 240; thus 1 + 2 + 3 + 4 + 5 -f 6 -f- 8 4- 10 4- 12 + 15 
+ 2() + 24 4- 30 + 40 f 60 - 240. 
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Perfect Numbers 



Constance Reid 

The Greeks, sreatly intrigued by the fact that the number 6 is the 
sum of all its divisors except itself (1 -f 2 -I- 3), called it a "perfect" 
number. They wondered ho%v many other such numbers there were. 
It was easy enough to ascertain by trial that the second perfect number 
was 28 (I + 2 -f 4 + 7 + 14). The great Euclid was able to prove that 
in all cases where a iunnf>cr can lie factored into the fomi 2" ' (2" — I) 
and 2 " — I is a prime numljer, the number must be the sum of all its 
divisors except itJH-'lf . Thus in the case of 6, n is 2 and 2" — 1 = S, a prime 
number: in the case of 28, n is 3 and 2" ~ I = 7, again a prime number. 
With Euclid's formula it was no difficult matter to compute that the 
third and fourth perfect numbers were 496 (n = 5) and 8,128 (n = 7). 
Hut l)cyond that the computation became laborious, and in any event 
it was not proved that this rule included all the perfect numbers. Euclid 
left for future mathematicians a challenging question: How many per- 
fect numbers are there? 

In more than 2.00(1 years mathematicians were able to turn up only 
12 numl^rs that met the strict re(|uirements for numerical perfection. 
Within the past year, however, the University of California mathema- 
tician R. M. Robinson has. with the aid of a modern computer, dis- 
covered five more. The discovery did not attract the attention of the 
press. Perfect nuinlK'rs are not useful in the construction of atomic 
bombs. In fact, they are not useful at all. They arc merely interesting, 
and their story is an interesting one. 

For many centuries philosophers were more concerned with the 
ethical or religious significance of perfect numbers than with their 
mathematics. The Romans attached the number 6 to Venus, because it 
is the product of the two sexes — the fxld (masculine) numla^ 3 and the 
even (feminine) number 2. The ancient Hebrews explained that Ciod 
chose to create the world in six days rather than in one because 6 is 
the nujre perfect numlK.'r. The eighth-century English theologian Alcuin 

(xiinted out that the second origin of the human race, from the eight 
nnnan licings on Noah's .\rk. was less |>erfcct than the first, 8 being 
an imperfect numl>er. In the 12th century Rabbi Josef Ankin recom- 
mended the study of |K.rfect ninnlH'rs in a program for the "healing of 
souls!' 

The mathematicians, meanwhile, had lx?cn making slow progress. 
The first four perfect numbers — 6. 28, 496 and 8,128 — had l>een known 
as early as the first centur)'. Not until 14 centuries later was the fifth 
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discovered. It was JiJI,550,JiS6 (n = 13). Then in 1644 the French mathe- 
matician Marin Merscnne, a colleague of Descartes, announced six more 
at one clip, and thereby litiked his name forever with perfect numbers. 
1 he numbers were now so lar^e that they were necessarily described 
only by the prime 4iumber 2" - 1. or, more bricHy. by the exponent, n. 
in Euclid's fonnula. The values of n for the 1 1 perfect numbers, includ- 
ing Mersennc's six new ones, were 2, 3, 5. 7. 13. 17, 19, 31. 67, 127 and 
257. In other words, the largest prime in the series was the enornmus 
number 2"' - 1. 

It was obvious to other mathematicians that Mersenne could not have 
tested for primality all the numbers he had announced. Kut neither 
could they. At that time the only method of testing was to try every 
p<jssible divisor of each number. By this laborious method mathemati- 
cians did test Mersenne s i'?rst eight numliers and found them prime. 

It was the great Swiss mathematician Leonhard Euler who tested the 
eighth number (2" - 1). Euler also proved that all even perfect numbers 
must Ik.' fif the form expressed by Euclid's theorem. No odd perfect 
number has ever been found, but it has never been proved that such a 
number cannot exist. 

For mi>re than 100 years the perfect numl>er formed from the prime 
2" - I remained the largest proved. Then in 1876 the French mathe- 
matician Eduard Lucas worked out a method by which a possible prime 
could Ik? tested without trying all potential divisors. At the same time 
he announced that he had tested 2'-' - I by his method and found it 
prime. 

According to Lucas, the numlK?r 2' I is prime if, and only if. it 
divides the (n 1 ) term of a certain series. In this series the first numlKr 
is 4 and each succeeding mnnlK-r is the scjuare of the preceding one 
minus 2; in other words 4. 14. I'JI. ^7.6M. and so on. For example, to 
test the prime number 7 (2 ' 1 ). one divides 7 into 14; the n - 1 term in 
ihis case being the second numlxr in the series, since n is 3. Since 7 
divides evenly into 14. it is prime by Lucas' tesi. 

Obviously even Lucas' short cut method iKcomes rather unwieldy 
when, as in the case of 2' ' ~ 1. one must divide 1 7»». 1 4 1.183,460.461: - 
2.1I.73l.687..303.715.884.nK'i.727 into the 12fith term of Lucas' series. 
For such lunnlKTS, mathematicians use a short cut of the short-cut: in- 
stead (A s«{uaring each term of the series, they srjuare only the remainder 
after they have divided the numln'r Inking tested into it. 

Even with the help of Lucas' method mathematicians were not able 
to finish testing all of the pissiblc Mersenne numl»ers until a few years 
ago. I heir tally showed that Mcrsenne's list of jx;rfect numljers was 
incorrect. He was right on nine innnlKrs (those for which n is 2. 3. f*. 
7. 13. 17. 19. 31 and j27). but he was wrong on two he had listed (those 
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with the exponents 67 and 257), and he had missed three numbers in the 
scries (with exponents 61, 89 and 107). Thus the list stood at 12, with 
2"'*(2"' — 1) the largest known i^rfect number. 

2(22 1) 
22(23 - 1) 
24(25 _ 1) 
26(27 - 1) 

212(213 _ 1) 

216(2" - 1) 
218(219 - 1) 
230(231 „ 1) 

260(261 ~ 1) 
288(289 _ 1) 

2ic^(2i07 — 1) 

2126(2127 _ I) 

2520(2521 ». 1) 

2606(2607 _ i) 
21278(21279 _ 1) 
22202(22203 _ J) 
22280(22281 _ 1) 

lAS'V j»f j)crl{< t numbers stands at 17. The lust five were mldeil by SVVAC. 

Then on January 30 last year Robinson fed the problem to the 
National Bureau of Standards' Western Automatic Computer, known 
briefly as SVVAC. This is a high-speed machine: it can do an addition 
of 36 binary digits in 64 millionths of a second. Robinsons' job was to 
break down the Lucas method into a program of the 13 kinds of com- 
mands to which the SW AC rcs|)onds. The job was complicated by the 
fact that, while the machine is built to handle numbers up to only 36 
binary digits, the numbers he was %vorking with ran to 2,300 such digits. 
It %v'as, he found, very much like explaining to a human being how to 
multiply IflO-digit numbers on a desk calculator built to handle 10. To 
tell SWAC how to test a possible prime by the Lucas method. 184 sepa- 
rate commands were necessary. The same program of commands, how- 
ever, could Ik.' used for testing any numlier of the Mersenne type from 

The program of commands, coded and punched on paper tape, was 
placed in the machine's "memory!" All that %vas then necessary to test 
the primal ity of any Mersenne number was to insert the exponent of 
the new numl>er as it was to be tested. I'he machine could do the rest, 
even to typing out the result of the test — continuous zeros if the numl>cr 
was a prime. 

37 



4] 



The first number to be tested was 2''* — 1, the largest of the 1 1 num- 
bers announced by Mersenne. Twenty years before it had t^n found 
not prime by D. H. Lehmer, who worked two hours a day for a year 
with a desk calculator to do the test. It happened that this evening 
Lehmer himself, now the director of research at the Bureau of Standards' 
Institute for Numerical Analysis on the IJ,C.L.A. campus, was in the 
room. He saw the machine do in 48 seconds what had taken him an 
arduous 70U and some hours. But the machine got exactly the same 
result, 

SWAC then continued on a list of larger pc •ii? -limes. Men^nne 
had said that all eternity would not suffice to test w nether a given num- 
ber of 15 or 20 digits was prime. But within a few hours SWAC tested 
42 numbers, the smallest of which had more than 80 digits. One by one 
it determined that they were not prime. Finally at 10 p.m. a string of 
zeros came up: the machine had found a new^ perfect number. Its prime 
was 2'** — K Just before midnight, 13 more numbers later, another 
prime came up: 2*"' — 1. In the decimal system this is a number of 183 
digits. 

The machine continued testing numbers when opportunity afforded 
during the next few months. Last June the number 2' ~ 1 was found 
to be prime. In Octol^r, concluding the program, it established as prime 
the numbers 2^^* — 1 and 2'-''^ — 1. The latter is the largest prime num- 
Urr, of any form, now known. 

I he perfect numbers of which these primes arc components are, of 
course, much larger— so large that in comparison with them conven- 
tionally "astronomical" numbers seem microscopic. Yet. by a proof as 
old as Euclid, mathematicians know that these numbers are the sum 
of all their divisors except themselves — just as surely as they know that 
6=1+2 + 3. 

I hey still do not know, however, how many perfect numbers there are. 
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FOR FUR I HER READING AND STUDY 

There is a vast amount «f literature which deals with all aspects of 
the theory of numbers. The references below are a few of the readily 
available sources for additional information on prime numbers and 
perfect numbers. 

PRIMES and COMPOSI FES; FACTORIZATION 

Bam., W. W. R. Mathftnatiml Rtrrfations and Essays. Marmillun, 1912; pp. 59-75. 
Barnf.tt, I. A. Somi' Ideas About \'umbt r Theory. National Council o{ Icarhen> 

of Mathematics, 1961 : pp. 4-1 7. 
BowMtu. H. AM> IJ(m>:R.s, |. K. Aritlnnctintl Exrursitms. Dover Puhlitations, 1961; 

pp. 106-113. 

Carnahan, VVai.tfr. Methotls for Systemaiiially Seeking Factors of Numbers. 

School Scirnrr and Mathematics, vol. 52. pp. 429-485 ^952). 
Carnahan, Walter. Prime Numbers in Sequences. S<hool Scirtue and Mathr- 

matics.sol 54. pp. 81.H-315 (1954). 
DAvrNPORT. H. The Higher Arithmrtir. Harper & Bros., 1960; pp. 9-39. 
C;rant. Haroi.d. The Prime Nunilier Theorem. Scripta Mnthematiea, vol. 20. pp. 

2.35-236(1954). 

f IJZUK, D. AND Ti cHMAN, Z. Elementary Bounds for the Number of Primes. Scripta 

Mathemativn, vol. 1 1. pp. 179-182 (1945). 
I.KiiMFR D. N. History of the Problem of Separatinf- a Number into Its Prmie 

Factors. Scientific Monthly, vol. 7. pp. 227-234 (1918). 
Lmiimer. D. N. Hunting Big C;ame in the I heoiy of Numljers. Hcriptn Mathc- 

mati(a,va\. 1. pp. 229-235 (1933). 
MAR.SHAi.1.. W. I.. Some Pro|)eities of Prime Numbers. The Peuta^mi, vol. 8, pp. 

5-8(1948). 

Ore, Ovstein. S'umher Theoiy and lis Histoiy. McGraw-Hill. 1918; pp. 50-85. 
Rauemacher, H. and loF.Piir/, (). The Enjoyvieut of Mathematics. Princeton 

University Press. 195: : Chapters 1,11. and 20. 
REiaiMANN, VV. J. The Eascination of Xumher.s. Mcthucn, 1957; pp. 46-63. 
.Steinman, I). B. .\ Setond Secpiel to Eratosthenes. Scripta Mnthematiea, vol. 22, 

pp. 79-80(19.56). 

PERFECT NUMBERS; MERSENNE NUMBERS 
.•\ri:hibaij>. R. C. Mersenne's Numlicrs. Scripta Mathematira, vol. 3, pp. 112-119 
(1935). 

Brow.n, Alan. Miihi)M'ife< t Numbers. Scriptn Miilhemntica, vol. 20, pp. 103-106 . 
(19.54). 

Dickson, L. E. Perfect and .\mitalile Numlx-rs. Sctentifir Monthly, vol. 10, pp. 

349-354 (April, 1921). 
Fraenkel, a. a. Perfect Numbers and ,\mi<able Numbers. Scnpla Mulhemalica, 

vol. 9, pp. 24.5-2.55 (1943). 
KRAVtTZ Sidney. Mcrsenne Numlx'rs. Recreational Mathematics Magazine, vol. 

8, pp. 22-24 (April. 1W2). 
McCarthy, Pai l. (Wd Perfect Numlxrs. Scripta Mathematua. vol. 23, pp. 43-47 

(1957). 
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Mersenme Primes. Robinson Primes, the I9th and 20th Perfect Numben. Recrea- 
tional Mathematics Magazine, \o\. 8. pp. 25-31 (April. 1902). 

Raoemacker, H. and Toeplitz. O. The Etijoymtnt of Mathematics Prinietun 
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Reio, CONSTANC4C. Fmm 7.ero to Infinity. London: Routledge k Keuan Paul, 1956- 
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loitJiARD. Jacxiijm. On Prime Numbers aiul Peifett Numbers. .SV»i>/fl Mathe- 

matica, vol. 1 9. pp. 35-39 ( 1 953). 
Uhlkr. Horace. A Brief History of the lnvestiJ^ation^ on Mersenne \uml>eis and 
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